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There have been long debates about the initial condition of inflationary perturbations. In this work 
we explicitly show the decay of excited states during inflation via interactions. For this purpose, we 
note that the folded shape non-Gaussianity can be interpreted as the decay of the non-Bunch-Davies 
initial condition. The one loop diagrams with non-Bunch-Davies propagators are calculated to 
uncover the decay of such excited states. The observed smallness of non-Gaussianity keeps the 
window open for probing inflationary initial conditions and trans-Planckian physics. 
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After decades of study of cosmic inflation, the initial 
condition for inflationary perturbations still remains a 
mystery. There is a long lasting debate on whether it is 
interesting to study non-trivial initial conditions, between 
the following points of view: 

• Inflation is an attractor solution. Thus all non¬ 
vacuum initial states should be washed away expo¬ 
nentially during inflation. So it is argued that only 
the lowest energy vacuum state is left, known as the 
Bunch-Davies (BD) vacuum [1]. 

• The non-vacuum states, once imposed, can stay for 
arbitrarily long time if the cosmological perturba¬ 
tions were linear. This is because, by definition, 
the linear perturbation theory relies on the initial 
condition in a trivial way. A modification of ini¬ 
tial condition implies different integration constants 
of the inflaton equation of motion, and thus gets 
directly imprinted on the cosmic microwave back¬ 
ground and the large scale structure [2]. 

There are indeed many reasons to take the second 
viewpoint seriously. The non-vacuum initial states for the 
inflationary perturbations are supported by a number of 
considerations: 

• Just enough inflation. If inflation starts not long 
before the largest observable scales exit the hori¬ 
zon, the initial condition of those largest observable 
scales cannot be the inflationary BD vacuum. Be¬ 
cause the BD vacuum is the lowest energy state of 
quasi-de Sitter background. Before inflation starts, 
the universe has a different background [3]. 

• Features during inflation. The features in the infla¬ 
ton Lagrangian, either one-off or periodic, kick the 
initial state of perturbations. The former excites 
all perturbation modes at an initial time, up to 
the energy scale defined by the sharpness of the fea¬ 
ture, and the latter continuously excite perturbation 
modes [1]. 

• Trans-A problem. During the de Sitter expan¬ 
sion, the physical wavelength of perturbation modes 


starts off at small length scales and gets stretched. 
If inflation lasts long enough, when tracking back in 
time, the physical wavelength of any perturbation 
mode starts shorter than the UV cutoff A of the 
IR effective held theory. As a result, new physics 
enters, such as new massive states or new opera¬ 
tors. The most famous problem of such kind is 
the trans-Planckian problem [6], because all known 
effective field theories have a universal cutoff Mp, 
above which quantum gravity effects kick in. 

• The selection of non-linear vacuum. The practical 
definition of the inflationary vacuum is the low¬ 
est energy state for inflationary initial conditions. 
However, ambiguity of this definition arises when 
this practical approach is applied to gravitational 
fluctuations at nonlinear order. This is because 
the definition of energy, and thus the lowest energy 
state, depends on the choice of time coordinate, and 
thus is gauge dependent. The best known effect of 
such kind in empty Minkowski space is the Unruh ra¬ 
diation, that accelerated observer finds him/herself 
in a thermal state. At non-linear orders for cosmo¬ 
logical perturbations, the curvature fluctuation C 
shifts the definition of time, and leaves ambiguity 
of the definition of vacuum in a similar way [5]. 

In view of the above considerations, the non-trivial 
initial states are considered to be an interesting possibility, 
and can be parameterized by the non-BD coefficients 
and C_(k), such that for any fluctuating field 0k 

0k = Mkflk + M^k^-k ) (1) 

where the mode functions are 

uk=c+(kvr+c_(k)(ur)*, (2) 

where A: = |k| and is the Bunch-Davies vacuum, prac¬ 
tically selected by the lowest energy state at early times. 
The consistency of uncertainty principle (commutation 
relation between 0k and its conjugate momentum) and 
the size of a quanta (commutation relation between Ok 
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and aj^,) requires 

|C+(k)|2-|C_(k)|2 = l . (3) 

Satisfying the above condition, the non-BD coefficients 
are then largely unconstrained for current researches, 
given that their back-reaction to the background energy 
density remains small [THS]- As a result, in the sub¬ 
horizon limit Iferj ^ 1, and in the weak non-BD limit 
|C'_(k)| <C 1, the inflationary two point function for the 
curvature perturbation C gets a correction 

(Ck(r)Ck'(T))' = (Ck(r)Ck'(T))'BD(l + (4) 


FIG. 1: One-loop diagrams are needed to illustrate the 
decay of non-BD coefficients. And such one-loop diagrams 
are related to the three point function squared through the 
optical theorem. 

seen in both of the two ways: decay products appear, and 
the non-BD state itself decreases towards the BD vacuum, 
as 


where prime denotes that (27r)^d^(k-|-k') is stripped, and 

Anon-BD('r) = -2Ck COs(2fcT -h 6»k) , (5) 

where C_ (k) is parameterized with real parameters Ck 
and dk as 

C_(k) = cke*®‘‘ . (6) 

The arbitrariness of choosing those non-BD coefficients 
seems to contradict the intuition that inflation is an attrac¬ 
tor solution and all non-ground states should be diluted. 

In this letter, we clarify the relation between the possi¬ 
bility of non-BD initial condition and the attractor feature 
of inflation. The key observation is that, the decay of the 
non-BD states, if possible, must happen through interac¬ 
tions. The process is similar to the process towards ther¬ 
mal equilibrium, where short modes radiate long modes 
and the system finally stabilizes at the equilibrium state. 

We start from commenting on the non-Gaussianities 
from the non-BD states. It is well known that the non-BD 
non-Gaussianities peaks at the folded limit [3 [TO]- For 
n-point correlation function, the limit is ki = ^^^2 
(starting from 4-point correlation function and for large 
C -, other folded configurations such as fci + fc 2 = ^3 + ki 
also have poles). The physical interpretation is simple 
(but not noticed in the literature to the best of our knowl¬ 
edge). Once a perturbation mode is initially not in its 
lowest energy state, it is unstable and can decay into 
longer wavelength modes. Thus the folded limit of non- 
Gaussianity simply captures the decay of those non-BD 
coefficients. 

To further clarify the situation, and make qualitative 
predictions, we would like to study explicitly how those 
nonlinear effects provide corrections to C'+(k) and (^-(k). 
For this purpose, we have to integrate over the decay 
products which we are not interested in, and consider 
the two point correlation function of the remaining high 
energy mode. This is nothing but the one loop diagram 
for the two point function, with non-BD coefficients. The 
relation between the non-Gaussianity and loop diagram 
is illustrated in Fig. [T] 

Formally, the above argument is just the optical theo¬ 
rem, as a requirement of unitarity in QFT - from proba¬ 
bility conservation, the decay of non-BD states must be 


cf = Ck exp [-r(T - To)] , (7) 

where tq is the initial time when the non-BD coefficients 
are generated, and F is the decay rate (which may depend 
on k though not explicitly written here). F can be time 
dependent, but the time dependence should be milder 
than linear on r — tq . 

In the remainder of the paper, we thus calculate the 
one loop correction to non-BD coefficients. To show the 
effect clearly, we choose the simplest possible model. Nev¬ 
ertheless, we expect that the decay of non-BD effect is 
completely general, and the method applies straightfor¬ 
wardly to other interactions. 

The model we choose is general single field inflation, 
where C = y/—gP{X, tp). Further we require that the 
sound speed is still c^ = 1, which is technically natural 
El from the effective held theory point of view, and 
provides us a simpler model for calculation. Namely, we 
require 

P,xx = 0 , A = X'^Pxx + ^X^P,xxx ^ 0. (8) 

After those simplihcations, the free Lagrangian for the 
curvature perturbation is 

£2 = ea^C^ - ea(5iC)^ , (9) 

and the 3rd order Lagrangian and Hamiltonian are 

c^ = -2a^^C\ n3 = -C3. ( 10 ) 

Following the standard QFT quantisation procedure, 
the C field can be written as 


Ck = UkOk + u 


-k“-k ’ 


( 11 ) 


with the mode function 

C+ (k) (1 + ifcr)e-*'=^ FC. {k){l -ikT)P'^^ . 

( 12 ) 

One can then calculate the one loop correction with 
non-BD coefficients making use of the in-in formalism. 
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With two interaction vertices, as illustrated in the left 
panel of Fig. the corrections can be calculated as 

(Ck(T)Ck'(r))3 = / dTi f (ir2(i//(Ti)Ck(T)Ck'('r)i?7(T2)) 

-'To J Tn 


-2Re / dTi / dT 2 {(w{T)Ck'{T)Hj{Ti)Hl{T 2 )) , 


(13) 


where the subscript 3 denotes that we are calculating the 
loop from C 3 . The contribution to equation (13) breaks 
into three parts: 

• Part without any non-BD coefficients C-, which is 
irrelevant for our purpose. 


contribution turns out to be a correction proportional 
to cos(0k + 2A:t), which corrects the amplitude Ck of the 
non-BD coefficient. On the other hand, the p + q — k > A 
part, including the UV divergence (which should be renor¬ 
malized), is proportional to sin(0k + 2A:r), which indicates 
a shift of phase 0k, and has no contribution to Ck at the 
leading order of the small Ck expansion. We are mainly 
interested in the decay of the non-BD amplitude, thus 
shall currently focus on the p + q — k < A part of the 
integral. We shall leave the full one-loop renormalized 
result in a future work |12] . 

After noting the above comments, the contribution near 
the apparent pole is straightforward to calculate. The 
leading contribution to the one loop two point function is 


• Part proportional to C_(q) or (^-(p), where q is 
the momentum running in the loop, and p = k — q. 
One should be able to find in this part that, once 
the loop contains a non-BD mode, how its energy 
is transferred to the other modes. This part of 
contribution does not directly show the decay of the 
non-BD mode. Mathematically, the contribution 
has an independent phase compared to 
Thus we shall drop this part at the moment and 
leave it to a future work m- 


• Part proportional to C- (k) . This shall be the part 
that we are going to calculate. 

Further, we take the sub-horizon limit |fcr| ^ 1. This is 
because at the super-horizon limit |fcr| 1, the pertur¬ 
bation is conserved and can no longer decay. The decay 
is dominated by the highest power of \kT\ terms, which 
are well captured in the |A:t| ^ 1 limit. 

Before showing the result, we would like to mention a 
subtlety here. Apparently, a divergence is encountered 
when we pick the highest order terms of r, by taking the 
T —> —00 limit. The momentum integral behaves as 


d^q 


p + q — k 


(14) 


Here the pole at p-\-q — k ^ Q arises because / dr exp[z(p-|- 
q—k)T] oc l/{p-\-q—k). As a result, equation (14) diverges 
logarithmly. On the one hand, the pole at l/{p + q — k) 
looks familiar, because it also appears at the calculation 
of non-Gaussianities, which shows evidence for the decay 
of non-BD. On the other hand, the contribution from this 
pole should not be physically infinity because otherwise 
it indicates an infinite decay rate, which is unphysical. 

This problem can be resolved by noting that near the 
pole, {p+q—k){T—To) may be finite even when kir—To) —>■ 
00 . Thus we should not consider (p + q — k){T — tq) to 
be always large in the |fcr| 3> 1 expansion. 

For this purpose, we split the momentum integral into 
p+q—k < A andp+q —fc > A parts, where A ~ l/(r—rg). 
The p + q—k < A part resolves the pole, where the contri¬ 
bution at p+g — fc is no longer singular. This p+q—k < A 


Anon-BD = A 


non—BD 


+ A 


1—loop 
non—BD 


(15) 


where 


A 1 —loop 
^non-BD 
A tree 
^non-BD 


cf (r) - Ck 


1 /A 


Ck 


3a2 


= -T2o UI 


1^ (^) +TTo+T^ f , 


(16) 


where we have chosen A « l/(r — rg) although there may 
be a coefficient discrepancy. 

It is interesting to note that this contribution is negative 
definite, indicating a decay of the non-BD coefficient. The 
full renormalization procedure may contribute to a change 


of coefficient in equation (16), but shall not change the 


sign or magnitude of the contribution. 

When /c(t — tq ) is very large, the one-loop contribution 
can dominate over the tree level contribution and the 
result becomes non-perturbative. And this case is indeed 
physically interesting because we are interested to study 
the limit where the non-BD coefficient significantly de¬ 
cays. Fortunately, one can resum the one loop diagram to 
address the non-perturbative decay of the non-BD coeffi¬ 
cients, by dividing t — tq into smaller intervals. Once the 
intervals are small enough, the contribution is guaranteed 


to be small because equation (16) is proportional to the 


length of the interval. This process is formally nothing 
but the flow of renormalization group (RG). To do this. 


we rewrite equation (16) into m 


.eff 


dlogc'k 


(r) « - 


40 


Prk^T^dT 


(17) 


where the integration constant should be determined to 
be c^®(to) = Ck- As a result, in terms of equation ([^, the 
decay rate is 


F = 


3 


(^) + + . (18) 


In terms of the non-Gaussianity estimator ~ If S ’ 
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the conformal decay rate is 


r = 


19683 

20000 



(19) 

Large non-Gaussianities have been pursued for two 
decades and a lot of people have been disappointed by 
the lack of large non-Gaussianity from recent experiments. 
However, here we show that large non-Gaussianity indi¬ 
cates faster decay of the non-trivial initial state of in¬ 
flationary perturbations. This is the price to pay for 
large non-Gaussianity. Nature is equally kind to keep 
non-Gaussianity small, which leaves a broader window 
open for probing non-standard vacua, and thus gives us 
more hope for probing how inflation gets started, features 
during inflation, trans-Planckian effects, and so on. 

For example, for I/^lI ~ 100, which is of order the 
current observational bound, the non-BD fc-modes which 
are 2 or more e-folds inside the horizon decay away signif¬ 
icantly. For I/nlI 1; which is the minimal requirement 
for the inflationary background dynamics being nonlin¬ 
ear, significant decay occurs at 4 or more e-folds. For 
I/nlI 0 . 01 , which is the slow roll bound given the 
detection of ^ 1 , the decay happens at 6 or more 
e-folds. The non-BD physics for modes with larger k are 
exponentially harder to probe. 

For the purpose of comparison, we also calculated the 
one loop diagram from the four point contact interaction, 
with Lagrangian 


A = (20) 

where /i can be obtained from the P{X, (p) model as 

M = ^X^P^xx + ‘^X^P^xxx + -^X^P^xxxx ■ (21) 

The diagram is illustrated in Fig. For this diagram, 
there is no fc-dependent non-BD pole as the £3 case. This 
is because, physically, this diagram is not related to the 
three point function by the optical theorem. Mathemati¬ 
cally, the external momentum k does not enter the loop. 


FIG. 2: The 4-point contact interaction does not contribute 
to the decay of non-BD coefficients at one loop level. 

To summarize, we calculate the one loop correction to 
inflationary two point function, with non-BD coefficients. 
This calculation uncovers the decay of non-trivial initial 
states for inflationary perturbations. The conformal de¬ 
cay rate of the non-BD coefficients are proportional to 


Thus smaller /nl indicates slower decay of 
the non-BD coefficients, and better preserves information 
from either beginning of inflation, features during infla¬ 
tion, or evidences of higher energy scales during inflation. 

The calculation can inspire a lot of future work. Here 
we only consider the dominant contribution from the 
simplest model of interactions. A full calculation of the 
current model is in progress | 12 j . and investigation of 
more general models are valuable. Also, it is interesting 
to go beyond the \C-\ <C |C+| assumption, which is 
currently imposed in our calculation. The non-BD mode 
may also decay into more than two softer particles, which 
corresponds to two or more loop corrections to the two 
point function of the curvature fluctuation. We also 
plan to investigate explicitly how the folded limit non- 
Gaussianity get regularized because of the decay of the 
non-BD vacuum m- 


Acknowledgments. We thank Xingang Chen for his many 
valuable suggestions at the early stage of this work. We 
thank Gary Shiu for helpful discussions. YW and HJ 
are supported by a startup grant from the Hong Kong 
University of Science and Technology. 


* Electronic address; phyw@ust.hk 

[1] T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. Bond. 
A 360, 117 (1978). 

[2] R. Easther, B. R. Greene, W. H. Kinney and G. Shiu, 
Phys. Rev. D 64 , 103502 (2001) hep-th/0104102 . 

[3] X. Chen and Y. Wang, JCAP 1407, 004 (2014) 
arXiv: 1306.0609 [hep-th]]. 

[4] X. Chen, JCAP 1012 , 003 (2010) arXiv: 1008.2485 [hep- 
th]]. 

[5] X. Chen and Y. Wang, arXiv: 1207.0100 [hep-th]. 

[6] J. Martin and R. H. Brandenberger, Phys. Rev. D 63 , 
123501 (2001) hep-th/0005209 . 

[7] X. Chen, M. x. Huang, S. Kachru and G. Shiu, JCAP 
0701 , 002 (2007) hep-th/0605045 . 

[8] R. Flauger, D. Green and R. A. Porto, JCAP 1308 , 032 
(2013) I arXiv: 1303.1430 [hep-th]]. 

[9] A. Aravind, D. Lorshbough and S. Paban, JHEP 1307 , 
076 (2013) I arXiv: 1303.1440 [hep-th]]. 

[10] X. Chen, B. Hu, M. x. Huang, G. Shiu and Y. Wang, 
JCAP 0908 , 008 (2009) arXiv:0905.3494 [astro-ph.CO]]. 

[11] D. Baumann and D. Green, JCAP 1109, 014 (2011) 
arXiv:1102.5343 [hep-th]]. 

[12] H. Jiang, et. al., work in progress. 

[13] H. Jiang, et. al., work in progress. 

[14] Note that we still need k{r — tq) > 1. Otherwise we are 
forcing interactions to take place at coincident time, and 
the other contributions will be missing. 








